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Abstract 



The COMPASS collaboration at CERN has investigated the ^ Tt^Tt^Tt^ reaction at center- 
of-momentum energy below five pion masses, y/s < 5m„, embedded in the Primakoff reaction of 
190 GeV pions impinging on a lead target. Exchange of quasi-real photons is selected by isolating 
the sharp Coulomb peak observed at smallest momentum transfers, t' < 0.001 GeV^/c^. Using 
partial-wave analysis techniques, the scattering intensity of Coulomb production described in terms 
of chiral dynamics and its dependence on the 37r-invariant mass ms^i = ^/s were extracted. The 
absolute cross section was determined in seven bins of y(s with an overall precision of 20%. At 
leading order, the result is found to be in good agreement with the prediction of chiral perturbation 
theory over the whole energy range investigated. 
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For hadronic interactions at low energy, quantum chromodynamics (QCD) can be formulated in terms 
of an effective field theory resulting from the systematic treatment of chiral symmetry and its breaking 
pattern, called chiral perturbation theory (ChPT). In this approach, the pions (7r+, 71°, tt^) are identified 
as the Goldstone bosons associated with the spontaneous breaking of the chiral symmetry. Properties and 
interactions of the pions provide the most rigorous test of ChPT as the correct low-energy representation 
of QCD. 

ChPT has accurately predicted the dynamics of low-energy nn scattering as measured in various kaon 
decay experiments, see Ref s. in |[T| . High-precision calculations of the pion polarizability in the ny — )• 
ny reaction have been performed, which however are at large variance with the existing experimental 
determinations Q. 

In view of this situation it is of great interest to study other low-energy reactions involving pions and pho- 
tons, apart from measuring the pion polarizabilities with improved accuracy as proposed at COMPASS- 
n |[3|. The ny — )• 371 reaction examined here allows for testing the chiral dynamics of pions in the region 
near threshold where calculations are most reUable ||T|. 

COMPASS at the CERN Super Proton Synchrotron is a large-acceptance, high-precision spectrome- 
ter Q, which is particularly well suited for investigations of high-energy reactions of particles at low to 
intermediate momentum transfer to fixed targets. The apparatus acceptance fully covers the phase space 
for reaction products emerging in forward direction. At very low momentum transfer in the process 
7r^Pb— X^Pb, photon exchange becomes important and competes with strong interaction processes like 
diffractive excitation via Pomeron exchange. The separation of these two contributions is an important 
aspect of this Letter. 

Primakoff reactions are peripheral hadronic reactions in the quasi-real photon field surrounding a heavy 
nucleus ||5J. The 7r-nucleus cross section can be connected to the 717 cross section using the well-known 
equivalent photon approximation (EPA) [61 



Here, the cross section for the process 7r^Pb— )• X Pb is factorized into the quasi-real photon density 
provided by the nucleus of charge Z, and Oy, the cross section for the embedded %^y ^ X scattering of 
a pion and a real photon. The mass of the charged pion is denoted m-j^, and JO,, is the ?i-particle phase- 
space element of the final-state system X . From the 4-momentum transfer squared t = {pji — pxY, the 
positive quantity = |f| — |f|niin is derived with |f|min = (■^ ~ '^l)^/(4£'bMm) for ^ given final-state mass 
mx = ^/s, where s = {pjt + py)^ is the squared centre-of-momentum energy. The lead form factor is 
approximated by the equivalent sharp radius formula F^f (?') = j'l {r^/f) with r = 6.84 fm. 

For the specific reaction studied in this Letter with X — 71^71^71+, the Feynman graphs for o^^"^ are 
depicted in Fig. [ija). The Dalitz plot distributions given in Fig. [TJb), which represent a 2-dimensional 
projection of the 5-dimensional phase space O3, show that the momenta of the 3 outgoing pions are far 
from being uniformly distributed in the m'^ji region near threshold. 

The data presented in the following were recorded in the year 2004 with a 190 GeV hadron beam (at 
the target composed of 96.8% 71^, 2.4% K , 0.8% p) interacting with a 3 mm lead target. The trigger 
selected events with at least two outgoing charged particles at scattering angles less than about 50 mrad. 
In the data analysis, exactly three outgoing charged particles, assumed to be pions, were required to 
form with the incoming beam particle a vertex that is consistent with an interaction within the target 
volume. Neglecting the tiny tai^get nucleus recoil at low t' , their energy sum £"3;!: equals the beam energy 
for the exclusive 7r^Pb— > 7r^7r^7r+Pb reaction. Taking into account the spread of the beam energy, 
Cbeara ~ 1-3 GcV, £"3;!: Is rcqulrcd to lie within ±4 GeV around the mean beam energy. Figure |2] shows 
the m-iji distribution with the requirement t' < 0.001 GeV^/c^, i.e. the Primakoff-f' region. The main 
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Fig. 1: (a) Leading-order, i.e. tree-level ChPT diagrams for n^y^ '^^ ^ embedded in the Primakoff reaction 
TT^Pb — > TT^TT^TT+Pb. (b) DaUtz plots corresponding to the leading-order ChPT cross section, normalized to the 
respective maximum of the distribution. The left plot is for ma^j- = 0.475 GeV/c^, the right one for 0.642 GeV/c^, 
i.e. at the lower and upper ends of the investigated region. The quantities and are the squared masses of 
the n^n'^ subsystems. 



contributions from diffractive excitation of the incoming pions to the ai (1260) and 7r2(1670) resonances 
are clearly seen. In the highlighted low-mass region, ms7t < = 0.72GeV/c^, leading-order ChPT is 
expected to be applicable [7]. The momentum transfer distribution for these low masses, excluding the 
kaon peak, is depicted in Fig. |3] The sharp increase with — )■ includes the Primakoff contribution that 
is investigated in detail in the following. 

The data are analyzed using an extension of the partial-wave analysis (PWA) package that was previously 
used for the high-f' events from the same data set [8]. TheX^ — n^n^n^ transition is modeled as a two- 
step process with an intermediate 2?! -resonance (isobar), X — )■ (2?! isobar)''?!^ — tt+tt^tt^, specified 
as 7^*'M'^ (271 isobar) [LJtt. Here J is the spin, P the parity and C the charge-conjugation parity of X, M 
the projection of J onto the beam axis, e the reflectivity of the final state in the Gottfried- Jackson frame 
(see e.g. ||8|), and L is the orbital angular momentum between the isobar and the bachelor pion. The 
partial waves found to be relevant at low m-^ji and low t' are summarized in Table [T] 

The new element of the present analysis is the replacement of the M= 1 isobaric waves at low masses by 
a single wave that is given by the fully differential form of the ChPT prediction and called from now on 
"chiral amplitude", see Eq. (17) in Ref. [1]. Replacing only M= 1 waves accommodates the transverse 
nature of the quasi-real photon as exchange particle, as M= is forbidden for real photons and strongly 
suppressed at small t' . This replacement is based on the observation [9] that Ti^Pb strong interaction is 
characterized by a (f')^exp(— Z7pi,f') behaviour at high energy and small t' , with the Pomeron slope pa- 
rameter ftpb ~ 400 (GeV/c)^^. Thus for strong Ti^Pb interactions at lowest t' , M= components dominate 
andM= 1 are suppressed, whereas photon exchange features a sharp spike at t' ?a \t\nim < 8- lO^^GeV^/c^ 
for masses mi,ji < m!,™, due to the t' dependence given by Eq. (fTl. 
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Fig. 2: Invariant mass spectrum of the final-state 3n system observed in the Primakoff region (f' < 
0.001 GeV^/c^), with the interaction point observed in the target region. The kaon decay into n^n^n^ is seen as 
sharp peak at niK ~ 0.493 GeV. The low-mass region of interest for this Letter is highlighted. 



Table 1: Wave sets used in the low-mass, low-f' region. In the figures showing PWA fits, in the low-mass region 
exclusively the "chiral amplitude" is used. In addition to the indicated waves, always a background contribu- 
tion that is homogeneous in phase space and the kaon decay contribution are allowed. Due to the very forward 
kinematics, for the M=l waves both reflectivities are allowed, as discussed in the text. 



M = 


M = \ 


0-+0+(7r7r)5 [5] TT 
0-+0+ p[P]n 
l++0+p[5]7r 
l++0+(;r7r)5[P]7r 
2-+0+(7r;r)5[D]7r 


1- +l^p[P]7l 

2++\^p[D]n 

2- +l^p[P]7i 
2-+\^{nn)s[D]n , 


f chiral 
[ amplitude 



In order to study the mass dependence of the partial waves in the Primakoff region t' < 0.001 GeVVc^, 
the low-mass (mjji < mg™) spectrum is subdivided into seven bins. In this region the PWA is performed 
with the isobar M= waves of Table 1 and the chiral amplitude for M= 1 . Above this limit both M= 
and M= 1 isobaric wave sets of Table 1 are used. In order to study the exchange mechanisms in this 
low-mass region, we look at the t' distribution extending the range up to 0.01 GeV^/c^. The intensities 
are obtained from a PWA with the M= set of Table 1 and the chiral amplitude for M= 1 both below and 
above 0.001 GeV^/c^. The results of these analyses are shown in Figs. |4|and|5| 

Figure|4]presents the fitted partial-wave intensities for the whole mass range, summed separately for M= 
and M= 1. The intensities are corrected for acceptance effects that are determined by a Monte Carlo 
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Fig. 3: Momentum transfer distribution of events in the investigated low-mass region. The Primakoff region 
(f' < 0.001 GeV^/c^) is highlighted. In order to suppress the kaon contribution, m^jc > 0.51 GeV/c^ is required. 



simulation and applied in the fitting process in their fully differential form in the 3-pion phase space. 
This correction does not deviate more than it 8% from the average value of 0.58 over the kinematic 
range considered. About 82% of the total intensity is attributed to M= waves, which are dominated by 
the <3i(1260) and 712(1670) resonances. 

For the low-mass region integrated between 0.51— 0.72 GeV/c^, Fig. [s] shows the PWA results for the 
t' distributions separately for the M=0 and ChPT (i.e. M=l) contributions. In both cases, an expo- 
nential function is fitted to the spectrum. For M= 0, Pomeron exchange with the slope bp^, as spec- 
ified above is confirmed. For M= 1, the exact shape expected for photon exchange can not be re- 
solved experimentally since it is distorted by multiple scattering and detector resolution. The result- 
ing spectrum is well described by an exp{—b't') function, where the theoretically "forbidden" for- 
ward region at f' — )• is filled by multiple scattering. The slope fitted to the ChPT intensity shown 
in Fig. [s] Zj' = 1447 it 196 (GeV/c)^^, is in good agreement with the simulation of pure photon ex- 
change, « 1600 (GeV/c)^^. Strong interaction M= 1 exchange, on the other hand, would have a 
slowly rising oc t' exp{—bt') distribution. From additional studies performed for the neighboring t' re- 
gion, 0.001 <t'< 0.01 GeV^/c^, strong interaction M= 1 contributions to the Primakoff region were 
estimated to be well below 5%. 

The fitted chiral amplitude very efficiently replaces all M= 1 waves, especially the dominant 1++1+ 
contributions also used as isobaric waves in PWA fit attempts of the low-mass region, see Table [T] 
Fitting the chiral amplitude as M= 1 contribution leads to smaller statistical uncertainties due to the 
smaller number of free parameters, while it features a signal strength compatible with that of the sum of 
the M- 1 waves in the isobaric description. The PWA attributes a comparable likelihood to both cases, 
which reflects that the real event distribution in the 5-dimensional phase space follows equivalently well 
the distributions expected from the two ansatze. But, the observed enhancement of M= 1 waves at small 
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Fig. 4: PWA intensities summed separately for the M=Q and M= 1 components. Total uncertainties are shown. 
Fluctuations of the systematic uncertainties, as discussed along with Fig. [6] were smoothed in the range 0.5 < 
rriijc/ (GeV/c^) < 1, and the obtained uncertainties were added quadratically to the statistical ones. 



t' lacks a theoretical explanation other than its chiral nature, and the comparison to the fit with isobaric 
waves is presented here only as a well-known starting point of the employed PWA. 

It is specific for the very forward kinematics studied here that the production plane spanned by the incom- 
ing beam and the outgoing X system becomes less accurately defined as t' — 0. As a consequence, the 
PWA not only shows the expected 1+ states, but also states with M^= I resulting from the limited 
resolution in the azimuthal orientation of the 37r-system, which applies only for M / 0. As partial waves 
are also identified by their signature in the remaining 4 kinematical variables that are not distorted by res- 
olution effects in the f' — )• limit, the two reflectivity contributions have been summed up and attributed 
to the well-motivated and expected M^= 1+ intensities. The validity of this approach was confirmed by a 
special Monte Carlo simulation, where the appearance of negative-reflectivity contributions arising from 
purely positive-reflectivity input states through resolution effects of the set-up was confirmed. 

In order to compare the observed M= I strength at mT,ji < m!,™ with the ChPT cross section prediction, the 
intensities obtained from PWA must be normalized to incoming beam flux and target thickness. While 
the latter is known to high precision, the effective beam flux is influenced by a variety of effects which are 
less well controlled. The beam intensity was not monitored with high precision, and absolute trigger and 
detection efficiencies have to be modeled. The method employed here to overcome these uncertainties is 
based on the free decay of the kaon component in the beam 1 10 1. As can be seen in Fig.|2j the observed 3- 
pion final state contains a contribution from kaon decays K — n^n^n^ . Since the fraction of kaons in 
the beam and their decay kinematics are known to a much better precision than the statistical uncertainty, 
the observed kaon strength was used to determine the effective luminosity. Uncertainties in the Monte 
Carlo description of the experimental setup cancel out, since they affect the K — Tt^n^n^ decay and 
the 71^7 n^n^n^ reaction in a similar way. 
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Fig. 5: PWA intensities for the investigated low-mass range. Each t' bin represents an independent fit to the data, 
yielding simulaneously the two values for the M= and the chiral intensity. 



In Fig. [6] the real-photon absolute cross section is shown. It has been obtained by taking into account 
the effective luminosity for the chiral amplitude and dividing out the quasi-real photon density factor 
as given in Eq. ([T]). A good agreement with the leading order ChPT prediction [1] is observed. This is 
expected since in ChPT loop effects and resonances are supposed to become important only at higher 
masses Q. 

The systematic uncertainties shown in the figures contain the following components: (i) the variation 
due to the fitting model; as different fit models gave results with similar values for the likelihood, we 
attributed the corresponding variations of the wave intensities as systematic uncertainty; (ii) the uncer- 
tainty originating from the luminosity determination, which is also indicated separately in Fig. [6} (iii) 
the full calculation of the e.m. radiative corrections is not yet available for n^j ^ n^n^n^, in contrast 



to the case of tt y ^ % n^n^. A first estimate indicates that this correction is well below 5% 1 11 1. All 
these systematic uncertainties sum up to about 18% in total. 

In conclusion we have measured the tt^/— )• n^n^n^ cross section in the region msji < Snijc, with a total 
uncertainty of about 20%. The measured points are in good agreement with the ChPT prediction. At this 
level of accuracy, ChPT therefore provides a good description of the strong pion-pion interactions at low 
energies, including the coupling to quasi-real photons via the Primakoff process. 

We gratefully acknowledge the support of the CERN management and staff as well as the skills and 
efforts of the technicians of the collaborating institutions. The authors thank N. Kaiser for instructive 
and helpful discussions. 
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Fig. 6: Cross section for n^y ^ k^k^k^ as a function of the total collision energy y(s = m^,i. The error bars 
show only the statistical error. The evolution of the systematic uncertainty, as discussed in the text and smoothed 
over the bins, is indicated as dashed line. 
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